Density modulation and electrostatic self-consistency in a two-dimensional electron 
gas subject to a periodic quantizing magnetic field 
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We calculate the single-particle states of a two-dimensional 
electron gas (2DEG) in a perpendicular quantizing mag- 
netic field, which is periodic in one direction of the elec- 
tron layer. We discuss the modulation of the electron den- 
sity in this system and compare it with that of a 2DEG in 
a periodic electrostatic potential. We take account of the 
induced potential within the Hartree approximation, and cal- 
culate self-consistently the density fluctuations and effective 
energy bands. The electrostatic effects on the spectrum de- 
pend strongly on the temperature and on the ratio between 
the cyclotron radius R c and the length scale 0$, of the density 
variations. We find that Os P can be equal to the modulation 
period a , but also much smaller. For R c ~ as p the spectrum 
in the vicinity of the chemical potential remains essentially 
the same as in the noninteracting system, while for R c <C a* 
it may be drastically changed by the Hartree potential: For 
noninteger filling factors the energy dispersion is reduced, like 
in the case of an electrostatic modulation, whereas for even- 
integer filling factors, on the contrary, the dispersion may be 
amplified. 



I. INTRODUCTION 

The interest in nonuniform magnetic fields, with spa- 
tial variations on a nanometer scale, has been stimu- 
lated by several recent experimental realizations, 
magnetic quantum wellsEI or magnetic superlatticesx 
In the quasi-classical regime of low magnetic fields, the 
theoretical investigations have concentrated, on the com- 
mensurability oscillations of the repiS|tiifity,Er1a which are 
equivalent to the Weiss oscillationslilaliil that occur in the 
presence of a periodic electrostatic potential. 

The quantum regime of nonuniform magnetic fields 
with a strong variation of the order of 1 T within a dis- 
tance of a fe.w hundred nanometers is now experimentally 
accessible.tj For this regime single-particle quantum me- 
chanical calculations, concerning the tunneling through 
magnetic barriers or the bound states in magnetic wells, 
have recentlv-heen performed by Peeters, Matulis, and 
Vasilopoulos.lljO Coulomb interaction effects have been 
discussed by Wu and UlloaEj who studied the electron 
density distribution and the collective excitations in a 
magnetic superlattice with a short period, comparable 
to the average magnetic length. They found that the 
periodic magnetic field gives rise to an electron-density 



modulation, which is reduced due to the counteracting 
induced electric field. 

In the present paper we consider a two-dimensional 
electron gas (2DEG) in a perpendicular magnetic field of 
the form B = Bq + B mo( i where Bq is a homogeneous 
part and B mo d is, in the plane of the 2DEG, periodic in 
one direction with zero average. We describe in detail 
the charge-density response to the periodic part and the 
effects of the associated electrostatic potential. It will be 
very instructive to compare this situation with the modu- 
lation by a unidirectional periodic electrostatic potential 
Vext(x) (we include the charge — e in the definition of the 
electrostatic potentials, which are therefore rather po- 
tential energies). The homogeneous part of the magnetic 
field is assumed to be strong enough so that a descrip- 
tion in terms of Landau levels is adequate; we denote by 
Iq = y/n/eBo and by uiq = eBo/m the magnetic length 
and the cyclotron frequency associated with the uniform 
field Bq. For a fixed mean electron density po, the num- 
ber of relevant Landau levels is of the order of the filling 
factor v = (2itIq)po and thus inversely proportional to 
the average magnetic field. 

We first note that classically a magnetic field - modu- 
lated or not - has no influence in thermodynamic equilib- 
rium because it drops out of the integral over momenta 
in the partition function (Bohr- van Leeuwen theorerrJiB). 
Especially a magnetic modulation does not lead to a 
position dependence of the electron energy, which re- 
mains mv 2 /2, and does not affect the equilibrium elec- 
tron density. In contrast one expects (e. g. from Thomas- 
Fermi-theory) that modulation by an electrostatic po- 
tential should lead also to a modulation of the density 
5p(x) = —(V ex t{x)/p)pQ 1 where p is the chemical poten- 
tial. In the classical limit, i. e. for low average magnetic 
field with Tiloq <C ksT, the density response of the 2DEG 
to electrostatic and to magnetic modulations is thus very 
different. A pure magnetic modulation does classically 
not give rise to electrostatic effects, whereas an external 
electrostatic potential is screened by the induced Hartree 
potential. 

However, in the quantum regime of low filling factors 
the two types of modulation affect the density in a very 
similar manner: Both lift the degeneracy of the Landau 
levels and lead to dispersive bands. The homogeneous 
part Bo of the magnetic field restricts the spatial ex- 
tent of the relevant wavefunctions to the order of the 
cyclotron radius R c — lo^/2np + 1 where np is the in- 
dex of the Landau level at the chemical potential (that 
is np is the largest integer smaller than vj g s \ we assume 
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spin degeneracy g s = 2 in this work). If rif is small 
the modulation is typically slowly varying on the scale 
R c and we can represent the Landau levels as functions 
e n {x), varying on the same length scale as the modula- 
tion. The width Ae n of these bands is of the order of 
TieAB/m or AV, where AB and AV stand for the am- 
plitude of the modulation of magnetic field and electro- 
static potential, respectively. In this situation we expect 
that the density is determined by a "local" filling fac- 
tor u(x) = g s f(e n (x)) 1 where /(e) denotes the Fermi 

n 

function. For temperatures satisfying fcgT <C Ae n , this 
must lead to a density modulation of order ((/ s /27t/q) for 
both electric and magnetic modulations. 

For an electric modulation it is known that, due to this 
strong effect of the Landau level dispersion on the density, 
the inclusion of the Coulomb interaction in the Hartree- 
approximation (which we will refer to as the electrostati- 
cally self-consistent system) changes drastically the spec^ 
trum of the system for low filling factors. Wulf et. ai.EZl 
found that, for filling factors not too close to an even- 
integer value, the self-consistent result corresponds to a 
nearly perfect screening of the modulating external po- 
tential and the Landau bands are flat within the order of 
kgT. If for even-integer filling factor the chemical poten- 
tial lies in a gap, the screening is much weaker although 
still considerable. For potentials strong enough to yield 
overlapping bands, |Ae„| > hu>o , the screening around 
even-integer filling factors becomes nonlinear, featuring 
two bands touching the Fermi level, and the width of the 
band np therefore locked to hio^. 

The width of the Landau band at the Fermi level, how- 
ever, is the basic—element for understanding transport 
measurements&u'Ea and a major aim of this work is to in- 
vestigate its behaviour for a magnetic modulation when 
electrostatic self-consistency is properly accounted for. 

The paper is organized as follows. In section |l| we 
describe our mod el a nd the self-consistency problem in 
detail. In section III we treat the case of low filling fac- 



tors, corresponding to a magnetic modulation varying 
slowly on the length scale R c . In section IV we discuss 



the regime of lower average magnetic fields, where the 
cyclotron radius is not small against the period of the 
modulation. The numerical results we present are ob- 
tained using the material parameters of GaAs, namely 
the effective mass m = 0.067m e and the dielectric con- 
stant k = 12.4. The average electron density is fixed to 
p = 2.4 • 10 11 cm -2 , chosen such that vB = 10 T, the 
period of the modulation is a — 800 nm and we consider 
values of Bq between 10 T and 0.1 T corresponding to 
cyclotron radii between 7 nm and 740 nm. 



II. DESCRIPTION OF THE MODEL 

We consider an idealized 2DEG confined to the plane 
{r=(x,y)} and subjecl to a magnetic field B(r) = 
(0, 0,B(x)) whichj3 Trli3 within the plane, is directed in 



z-direction, does not dependlll on y, and has a simple 
periodic dependence on x : 



B(x) = B + Bi cosKx , 



(1) 



where K = 27r/a is the wave vector of the modula- 
tion. We start with the noninteracting 2DEG described 
by the standard single-electron Hamiltonian H° = (p + 
eA) 2 /2m in which we use the Landau gauge for the vec- 
tor potential, A(x) = {0,B x+ (B 1 /K)sinKx,0). The 
eigenf unctions of H° depend on y only through a plane- 
wave pref actor, 



i> nXo (x,y) = L-^ 2 e- lX "y/ l ^ nXo (x), 



(2) 



with L y being a normalization length and Xq the center 
coordinate. The functions 4> n x i x ) are the eigenvectors 
of the one-dimensional Hamiltonian 



H°(X ) = huo 



2 dx 2 



1 

w 



X 



s \ 2 

— sin Kx 
K 



(3) 



where s = Bi/Bq will be refered to as the modulation 
strength. 

For the homogeneous system, s = 0, the functions 
4>nX a {x) are oscillator wave functions centered on Xq, 
also known as Landau wave functions <p%x > associated 
with the degenerate Landau levels = (n + l/2)7k^o- 
The degeneracy is lifted for s ^ 0, and the resulting en- 
ergy bands e n (Xo), together with the corresponding wave 
functions, can be obtained by diagonalizing the reduced 
Hamiltonian (^) in the basis of the Landau wave func- 
tions. The matrix elements can be written as 



(<p L nXa | H°(X ) | <(%, Xo ) = 

E nn '(z) + \frm' E n ^,n' 

-V(ra+l)(n' + l)#n+i,n'+i(z) cos (kX + ( n -n')|) 

+ £ l 5nn ' ~ S ""' (4z) C0S ( 2KX ° + } ' (4) 

where z = (KIq) 2 /2. We have used the notation: 

/,\ 1/2 



E nn ,(z) = (^pj e~>l 2 z^y 2 L r ; L 7 n \z) 



= (-!)""« E n , n (z) 



(5) 



with L™(z) being a Laguerre polynomial. Applying first- 
order perturbation theory we get from (Q) the energy 
levels as simple cosine-shaped bands 



(X ) = hu (n + 1/2 ) + s G n {z) cos(KA ) (6) 

where the factor 2shuj^G n (z) = sTwjq e^ z / 2 {2L] l {z) — 
L°(z)) has an oscillatory dependence on the ratio 
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la\/2n + l/a which is the basic reason for the commensu- 
rability oscillations seen in transport experiments. Ella The 
limits of validity of Eq. (^) will be discussed below. 
The single-particle density is given by the formula 



+00 



p(x) 



9s 



2nl 2 

u n=0 



dX f(e n (X )) I <j> 

uXq (^) 



(7) 



where /(e) denotes the Fermi function and g s — 2 ac- 
counts for spin degeneracy. 

The density determines the electrostatic (Hartree) po- 
tential, which we treat by Fourier expansion V H {x) = 
E V" cos(rjKx). Here 

»7>1 



V H — — — — — 



1 



Pn 



r] 2tt as 



(2-Kll p„) Tilu , (8) 



with p(x) = ^2 Pri cos(t]Kx) and as the effective Bohr 
v>o 

radius. For GaAs, 2-Kas ~ 63 nm. We assume that 
the system is electrically neutral such that the average 
density po does not contribute to V H but only determines 
the chemical potential contained in the Fermi function. 
The Hartree potential has to be added to the Hamiltonian 
ml) and gives a contribution 



(tix I V H (x) I ri, Xo ) 

, {ifz) cos (r)KX + (n - n') |) (9) 



to the matrices (Q). The strongest influence on the in- 
duced potential originates from the low Fourier com- 
ponents of the density, which are related to the long- 
range charge fluctuations. We diagonalize the Hamilto- 
nian H°(X n ) + V H self-consistently with Eq. (0) by a 
numerical iterative scheme. 

To understand the way the system achieves self- 
consistency, we occasionally consider also a 2DEG sub- 
ject to a homogeneous magnetic field B and a cosine 
electrostatic potential 



ViO) = Vi cos(Kx) 



(10) 



with modulation strength v\ = V\/Tiojq instead of the 
magnetically modulated system (|l|). First-order pertur- 
bation theory yields for the electric modulation ([H]) the 
spectrum 



PTl, 



(X ) = htoo (n + 1/2) + vtFniz) cos(KX ) 



(11) 



with F n (z) = e~ z ' 2 L n {z). 



III. THE LIMIT OF LONG PERIOD 

In this section we deal with a long-period magnetic 
modulation, with a strong average magnetic field, such 



that KIq <C 1 and R c ~ lo (but not necessarily with 
B\ <C Bo). Approximate analytical results for both elec- 
tric and magnetic modulations will be developed for a 
better understanding of the energy spectra and electron 
density. We first describe the properties of the noninter- 
acting system. 



A. Noninteracting electrons 

Treating the magnetic modulation as a perturbation 
presents some difficulties in this limit, since for nonvan- 
ishing s and z — * the matrix elements {<PnX I H°(Xo) I 
Vn+mXn) gi yen by Eq. (Q) diverge for m = 0,±1, while 
those with m = ±2 are finite and those with | m | > 2 
vanish. Thus the Hamiltonian matrix becomes band- 
diagonal, and the divergent elements cancel in a compli- 
cated way in order to yield finite eigenvalues. Therefore, 
for KIq « 1 an accurate numerical diagonalization re- 
quires a large matrix (|]) and the Landau level mixing 
is strong, except if s — > 0. This complication does not 
occur for the electric modulation (10) for which the Lan- 
dau wave functions diagonalize the matrix in the long- 
period limit for any v\ , and first-order perturbation thp=. 
ory gives the exact energy spectrum for z — » 0, namelycJ 
£nx = (n + l/2)huj + Vcos(KX ) . 

Instead of using standard perturbation theory with re- 
spect to the modulation stregth s, we can handle the 
Hamiltonian (^ by performing a Taylor expansion of the 

around 



potential term (jiujQ /2l^j (x — Xq 
its minimum X\ given by 



iKxY 



X 1 = X ~-smKX 1 
K 



(12) 



For fixed KXq and |s| < 1 this has a unique solution 
KX X with X x = X Q for KX = 0,tt. The parabolic 
approximation reads 



H°(X ) « hw 



2 dx 2 



(1 + scos KX-i) 
2Zg 



-{x 



X,f 
(13) 



with an error term of order s TiuiqKIq({x — Xi)/IqY from 
which we can show that (|l3|) yields the eigenvalues and 
the density for low filling factors correct to leading or- 
der in KIq. The Hamiltonian ( |l3|) is equivalent to the 
unperturbed one, Eq.(|^), but with modified center coor- 
dinate Xi, cyclotron frequency Cjq = loq(1 + scos(KXi)) 
and magnetic length l a = lo/yjl + s cos(KXi). The 
main effect of the mag netic modulation on the wave func- 
tions is the shift (n2h of their center of weight. We see 
that, since K{X\ — Xq) is independent of K, the abso- 
lute shift X\ — Xq increases with increasing modulation 
period (K — > 0) at equivalent positions within the pe- 
riod (i. e. for fixed KXq) except for KXq — 0, tt. This 
explains the difficulties with the standard perturbation 
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theory which expands the shifted Landau functions with 
center Xi(Xo) in the basis of Landau functions centered 
around Xq. 

The Landau bands resulting from Eq.(^3|) are 



»(lo) = Rwb(l + a cos (KXiiXo))) (n + ~ 



(14) 



The appearance of X\ instead of Xq in Eq. ( |14] ) leads to a 
substantial deviation of the simple cosine band shape pre- 
dicted by first-order perturbation theory; the bandwidth 
is, however, given correctly by Eq. ([|). In calculating 
the density, the indicated replacement of Iq by Iq leads to 
corrections of higher order in Klg and, since this order 
is not included correctly, is not to be used. We there- 
fore insert just shifted Landau-functions into Eq. (ffl) and 
obtain 



p(x) 



^2 E/ dX f(e n (X )) | ti fXl{Xo) {x) 



2ttZ, 
;2' E 



dXt- 



dXn 



2nl 2 ^ J dXi 



/(e„(X (*i)) \rt tXl (x) | 2 . (15) 



From ( |12[ ) we have dX /dXi = 1 + scos(KXi) and with 
( |l4| ) for the energy spectrum we finally get 

p(x) ^Y I pY,f dXl ( 1 + s C0S ( KX i)] 

x/(awb[l + scos(A-Xi)](n + l/2)) K^z) | 2 . (16) 

Both results (0) and @ turn out to be reliable within 
a relative accuracy of (s KR C ) . They can also be derived 
by a simple variational approach, using a set of translated 
oscillator states tp^ x +u( x ) as trial wave functions and 
taking the limit z -^j-Q after minimizing the expectation 
value of the energy.Ej The numerical results which we 
shall present are obtained from a diagonalization of (^), 
however. 

We note that for the electric modulation ( liol) t he re- 
sults corresponding to Eqs. (|l2]), ( pT| ) and (p.6[) read 



X\ — Xq 
e n {X ) 



viKl^sin KXi 



TlLJQ 



1/2 



vi(l - (1/2) (n+ l/2)( J ftT/ ) 2 ) cos(.KXi) 



(17) 



and 

p(z) 



I dXx^l-vxiKhfcosiKXx 

n 

/ (nc; [n + 1/2 + «a cos(M )] ) | (a:) | 



(18) 



The error term in the Taylor expansion around X\ is here 
of order v\ hcuo(Klo) 3 ((x—Xi)/lo) 3 and permits inclusion 
of the (KIq) 2 terms. The total width of the bands from 



( |l7| ) is also obtained from the result ( pH ) of perturbation 
theory by expansion around Kl Q = up to order (Kl ) 2 . 

We see from Eqs. ( |i"7| , |l8| ) that for a long-period cosine 
electric modulation the bands follow the potential with 
constant width and the states are not changed by the 
modulation; consequently the density is only affected by 
the dispersion of the levels via the argument of the Fermi 
function. In contrast, for a magnetic modulation accord- 
ing to ( pT[ ) the widths of the Landau bands increase lin- 
early with n and the ^-dependent prefactor of the Fermi 
function in ( |l6| ) does not decrease with increasing period. 

In Fig. |l| the dashed lines show the modulation of the 
density of the noninteracting system with a magnetic 
modulation of amplitude B± = 0.1 T for different values 
of the filling factor between 4 and 6 obtained by sweeping 
Bq ; the temperature is 1 K, so that ksT is much smaller 
than sTiluq. The density is given in units of 1/(2ttIq) so 
that the mean value of each line equals v. The lines for 
the even-integer values of the filling factor, v — 4, 6, are 
marked with circles. They show a cosine form, where 
the amplitude is larger for v = 6 than for v = 4. This 
behaviour is easily derived from Eq. (|l6|); since here the 
chemical potential lies in a gap, the Fermi function is 
either or 1, and the integral gives to leading order in 
Klo 



6p(x) 



(s^/27t/q) cos Kx 



(19) 



In the corresponding result for the electric modulation, 
the factor s is replaced by —vi(KIq) 2 which has a differ- 
ent sign and vanishes for Klo —> 0. The persistence of 
an finite density modulation at even-integer filling factors 
for a period much longer than the magnetic length consti- 
tutes a major difference between the two types of modula- 
tion for strong average magnetic fields. Note that the re- 
sult jl9| ) can also be written in the form p{x) — v/2itl 2 (x) 
(y small and even) where l{x) = yJ%/eB(x) is the mag- 
netic length corresponding to the local field B(x). This 
means that we can in the long-period limit think of the 
magnetic modulation as changing the local degeneracy of 
the Landau levels, thus leading to a modulated density 
even for spatially constant filling factor (in this work we 
use the notion of an x-dependent filling factor v{x) as just 
counting the number of locally occupied bands, which 
makes sense of course only in the long-period limit). 

The dashed lines between the ones with circles in Fig. |l| 
show the behaviour of the density while the n = 2 level is 
successively filled. Due to the energy dispersion (|lj) and 
the low temperature, the n = 2 states around KXq = tt 
are occupied first, forming a region with local filling fac- 
tor v(x) = 6 while around KXq — 0, 2ir we still have 
v{x) = 4 until the total filling closely approaches 6. Since 
the spatial extent of the wavefunctions is small compared 
to the period a, the difference in density between these 
two regions is of order <7 s /(27tZq). We observe, however, 
that within a region of constant local filling factor the 
density is not constant but follows the cosine shape im- 
posed by Eq. (fL9j) with v replaced by the appropriate local 
filling factor vyx). 
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B. Self-Consistent System 

Since the density profiles of the noninteracting system 
correspond, according to Eq. (@), to electrostatic poten- 
tials with amplitudes larger than hu>o, we expect substan- 
tial changes in the spectrum when we take electrostatic 
self-consistency properly into account, as we do now. The 
resulting densities are plotted as solid lines in Fig. |l| and 
show much smaller fluctuations. In Fig. |we show results 
for a magnetic modulation of B\ =0.1 T at filling factors 
(a) v = 5, (b) v = 4 and (c) v — 14.3 corresponding to 
average fields B = 2.0 T, 2.5 T and 0.7 T, respectively. 
The upper panel displays the self-consistent spectra for 
temperatures T = 1 K and T = 0.1 K together with the 
noninteracting spectra, the lower panel shows the cor- 
responding self-consistent densities. More data for the 
self-consistent bandwidths and the density amplitudes in 
this regime are also displayed in Fig. || and |^ which are 
discussed in section |IVB . 



When the total filling factor is small and not too close 
to an even-integer value, the regions of increased den- 
sity correspond to the minima of the Landau bands (cf. 
Fig. [l]). Therefore the generated Hartree potential, which 
is maximum at maximum electron density, will act to re- 
duce the dispersion of the not fully occupied band with 
index np. The self-consistent solution yields then a very 
flat ( "pinned" ) band with deviations of only the order of 
kgT from the chemical potential, and the local filling fac- 
tor is fractional over the whole period. This situation is 
shown in Fig. (a) for an odd-integer average filling fac- 
tor v — 5. The self-consistent potential here has to cancel 
the dispersion of the not fully occupied level np, which 
is larger than the dispersion of the levels with n < rip. 
Thus the potential is V H (x) « — s Twoq (np + 1/2) cosKx 
and the dispersion of the levels with n < np is reversed 
in sign. 

For even-integer filling factor, however, according to 
Eq. (|l9|) the regions of increased density correspond to 
maxima of the Landau bands, since both p(x) and e n (Xo) 
follow the shape of the magnetic modulation with posi- 
tive sign. Consequently, the potential generated by the 
density modulation ( |19[ ) increases the dispersion of the 
highest occupied band np instead of acting against the 
modulation broadening. If the magnetic modulation is 
sufficiently weak, the resulting self-consistent potential 
can be calculated by combining Eqs. (|l6|), ( [l8| ) and (||) 
as V H (x) — Vh cosKx (y small and even) where 



V H = 



l + w(Kl ) 2 v 



(20) 



and w = a/2naB 1. This linear behaviour breaks 
down, however, if the resulting bandwidth |Ae„ F | = 
2shuJo(np + 1/2) + 2Vh exceeds hui . In this case the 
next-higher band reaches the chemical potential around 
KX = 7r and the self-consistent solution (shown in Fig.|| 
(b) for v = 4) features a region around KXq — 0, 2tt 
where the band np is pinned to fi, a region around 



KXq = 7r where the band np + 1 is pinned to /i and 
a region in between where the chemical potential lies in 
a gap and the density still follows the cosine shape jlS|). 
As described in section Q, similar effects of electrostatic 
self-consistenpj£ are obtained for an electrically modu- 
lated systemjUl but there the bandwidth of the highest 
occupied band is always reduced compared to the non- 
interacting results. Then, a modulation strength v± > 1 
is needed to produce the formation of pinned regions at 
even-integer filling factors, whereas in the magnetic case 
only sw ~ 1 must be satisfied. 

For the parameters of Fig. || (c) the bands of the non- 
interacting system do overlap at the Fermi level due to 
the linear increase of their width with n. In this situa- 
tion the density fluctuations and the induced potentials 
consist mainly of higher Fourier components. The cor- 
responding wavelengths 2-K/nK, with n > 1, are compa- 
rable to R c , even though R c <C a is still satisfied. We 
therefore cannot discuss the effects of the Hartree poten- 
tial here within the limit of a long period but instead we 
have to consider the density response for electric modu- 
lations with a ~ R c . This is done in the next section. We 
observe, however, that in Fig. || (c) the spectrum around 
the Fermi level remains essentially unchanged although 
we can tell from the behaviour of the lowest level that a 
considerable electrostatic potential does exist. 

For the lower temperature T = 0.1 K the density traces 
m Fig. | (b) and very pronounced in (c) show also su- 
perimposed short-period oscillations. These have their 
origin in the nodes and maxima of the wavefunctions 
and can also be reproduced by Eq. ([!(]) with the Lan- 
dau functions fn( x )- 



IV. R c COMPARABLE WITH PERIOD 

In this section we discuss properties of the modulated 
system obtained when for fixed modulation amplitude 
Bi the average field Bq is lowered so that we enter the 
regime where R c is no longer small compared to the pe- 
riod a. In this case the approximation of the modulation 
by the first terms of a Taylor series breaks down and its 
actual functional form becomes important. However, the 
numerical method outlined in section [n] is still valid pro- 
vided that s < 1, i. e. the total magnetic field B(x) is 
nowhere vanishing. We consider first the noninteracting 
system. 



A. Noninteracting Electrons 

The quantity we are most interested in is the ampli- 
tude Ae np of the Landau level np at the Fermi energy. 
In Fig.|| (a) (solid line) this bandwidth is shown for a 
weak magnetic modulation B\ = 0.01 T and average 
fields B Q in the range 10 T > B > 0.125 T. Starting 
from high fields at Bq = 10 T we have first np = 



5 



and the bandwidth is sTlluo. When the field is lowered, 
|Ae„ F | increases in steps of 2shio at even-integer filling 
factors, that is when np jumps by one, as follows from 
(fj). When2i? c / a becomes larger than about 1/4 the in- 
crease of | Ae„ F | becomes visibly slower and goes over into 
an oscillatory behaviour with the first maximum at about 
2R C = 0.6 a. This can be understood with the result (||) 
of first-order perturbation theory. Using the asymptotic 
relation between Laguerre polynomials and Bessel func- 
tions we obtain lor the bandwidth at the Fermi level from 
(0) the formulaQ 



|Ae, 



2 s Tiloq A m Ji (KR C 



(21) 



where J\ is the Bessel function of order 1, R c = 
l Q \/2n F + 1, and A rn = {R c /Kll). The expression 1211) 
describes well the bandwidths (also for small filling fac- 
tors) as long as the average field is strong enough to en- 
sure s -C 1. It has zeros at approximately 



KR C = (A + 1/4)tt , A = 1,2,. 



(22) 



corresponding to a.-flat-band with vanishing dispersion 
at the Fermi level.El For our parameters we encounter 
only the first (A = 1) of these magnetic flat-band situ- 
ations around Bq 1 = 6.2 T _1 . From Eq.([2l]) we infer 
further that the maximum values of the bandwidth at 
the Fermi level are of order 2shujoA m rather than 2s?uvo- 
If we replace in R c the discrete 2np + 1 by v, the factor 
A m becomes y 'a 2 pQ/irg s and depends thus only on the 
period and mean density; typically we have A m > 10, 
e. g. for our parameters A m = 15.6. Consequently a 
seemingly weak modulation strength s ~ (1/A m ) <C 1 is 
sufficient to yield for KR C around a maximum of the 
Bessel function J\ a bandwidth |Ae„ F | > Hloq which 
means that the bands around the Fermi level do over- 
lap. In Fig. ^ (a) the bandwidth at the Fermi level for a 
modulation of B\ =0.1 T is plotted; it is larger than Tvujq 
for Bo < 1 T. For this stronger modulation the spectra 
show, at fields B < 0.25 T, also substantial deviations 
from the first-order perturbation expression (^J) , because 
the modulation strength s then becomes too large. The 
bands around the Fermi level are not cosine-shaped in 
this regime but have extrema away from KXq — 0, n. 
As a consequence, the bandwidth does not go through 
a zero at the flat-band condition (^2|) although its be- 
haviour still resembles the oscillations described by (21). 
For an electric modulation the result corresponding 



to 



© is |Ae, 



2viTujJqJq{K R c ) 



with zeros at 



KR C = (A - l/4)7r , A = 1,2,... (these are the elec- 
tric flat-band situations) and the bandwidth is always 
smaller than 2v\hu>Q. In Fig.(^) (a) the bandwidth at 
the Fermi level for an electric modulation of amplitude 
Vi = 0.27 meV = A m (Ke/m) ■ 0.01 T is shown as dashed 
line; the modulation amplitude is chosen such that the 
bandwidths are comparable to the ones induced by the 
magnetic modulation also depicted in this figure. 



In describing the induced density fluctuations Sp(x) — 
p(x) — po we face the difficulty that these have in general 
no simple shape (see Fig. [j]). As a measure of their mag- 
nitude we therefore concentrate on their amplitude |Ap|. 
This quantity is displayed in units of (1/2ttZo) for the non- 
interacting case at temperatures T = 1 K and T = 0.1 K 
in Fig. |3| (b) for a weak magnetic modulation not leading 
to band overlap and in Fig. ^ (b) for a stronger modu- 
lation with overlapping bands. The lower temperature 
T = 0.1 K corresponds to k B T = 8.6 • 10~ 3 meV which 
in the displayed range of Bo can be considered as small 
compared to Tiujq whereas for the higher temperature 
T = 1 K the finite size of k B T becomes important for 
about Bq 1 > 3 T _1 . The results for the density can be 
summarized as follows: If the bands around the Fermi 
level do not overlap and k B T is small against Hloq and 
also against the gap between the bands np and np + 1, 
then for even-integer filling factor the density is cosine- 
shaped Sp(x)\ v cvcn = sr m (KR c ) cos(Kx). This is due to 
the distortion of the occupied wavefunctions by the mod- 
ulated magnetic field. The amplitude r m is for KR C — > 
equal to the total density po (see Eq. (|l^)) and shows for 
lower fields oscillations in KR C with zeros at both elec- 
tric and magnetic flat-band situations. During the filling 
of each band, i. e. when the filling factor is not an even 
integer, an additional density fluctuation of the order of 
(1/2ttIq) is produced due to the dispersion of the bands, 
like in Fig. |l|. Both effects have a tendency to cancel 
each other. If k B T is not small compared to Hlvq or to 
the gap between the bands uf and n F + 1, this cancel- 
lation becomes almost perfect and the resulting density 
modulation is minute. If the bands around the Fermi 
level do overlap (as in Fig. [| (b) for Bq < 1 T), the den- 
sity has a complicated shape with several extrema and 
its amplitude shows an irregular dependence on the fill- 
ing factor. For k B T -C hu the amplitude of the density 
fluctuations is still of the order of 1/2itIq (but not larger) 
whereas for a higher temperature again no appreciable 
modulation of the density is produced. 

We see thus that the modulated magnetic field affects 
the density only if R c is small compared to the modula- 
tion period or if the temperature is very low; in any case 
the resulting density modulation is limited in amplitude 
by (7 s /27tZq. For not too low temperatures the Thomas- 
Fermi-prediction, namely no modulation of the density, 
holds to good accuracy as soon as R c > a/4. 

The density amplitudes resulting from an electric mod- 
ulation are shown in Fig. || (c). As discussed in section 
III A, we find for R c < a/8 density fluctutations of order 



g s /27r^Q if the level rip is partially occupied, and essen- 
tially no modulation of the density at even-integer filling 
factors. For lower fields Bq, however, the density modu- 
lation becomes dominated by a cosine contribution whose 
amplitude is not related to <7 s /27r/g but rather equals the 
Thomas-Fermi value 

Sp TF (x) = -(V(x)/p)p . (23) 

Since Sp TF is independent of Bq it appears in the plotted 
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quantity 27tZq|Ap| as a linearly increasing background. 
For the higher temperature we find that for R c > a/4 the 
density is described accurately by Eq. (p3|), i. e. Sp(x) = 
— (vihtuo/(J-)po cos(Kx) . For the lower temperature devi- 
ations from this result appear which are of similar mag- 
nitude as the corresponding deviations from zero for the 
magnetic modulation (except that at even-integer filling 
factors they do not vanish at the magnetic flat-band con- 
dition but only at the electric one). The density mod- 
ulation for the low temperature and even-integer filling 
factor is entirely due to the distortion of the occupied 
wavefunctions and follows, very well the formula derived 
by Aleiner and GlazmanEH from first-order perturbation 
theory. Most important for our purposes is the fact that 
in any case for R c > a/4 there exists an appreciable mod- 
ulation of the density, whose main part is a wave-function 
effect and follows the electrostatic potential linearly, in- 
dependent of the spectrum at the Fermi level. Due to 
the linearity in V, it is clear that this applies also to the 
density response to those higher Fourier components of a 
non-cosine electrostatic potential whose wavevectors r\K 
satisfy r\KR c > 1 . 



B. Self-Consistent System 

Having discussed the density response induced by mod- 
ulated magnetic and electric fields, we now proceed 
with the investigation of the effects of electrostatic self- 
consistency for the magnetically modulated system. In 
Fig. [| and ^ (a) and (b) the bandwidth at the Fermi level 
of the self-consistent system and the amplitude of the 
Hartree potential | AV H | are shown for the same param- 
eters as used in Fig. || and |[ respectively. The bandwidth 
of the noninteracting system and fuvo are also shown for 
comparison. We can clearly distinguish the high-field 
regime (limited by R c < a/4) discussed in section [II 
where the spectrum around the Fermi level is dominated 
by the electrostatic effects: Instead of the monotonous in- 
crease in the noninteracting system the bandwidth at the 
Fermi level is here of order ksT when the filling factor is 
not close to an even-integer value and has at even-integer 
filling factors sharp maxima with a height of the order 
hujQ (this value is reached only for the stronger modula- 
tion in Fig. H while in Fig. || Eq. ( |2fj| ) remains valid). The 
amplitude of the potential equals the difference of the 
self-consistent and noninteracting bandwidths, reflecting 
the fact that electrostatic self-consistency is achieved by 
adjusting the spectrum. 

For lower fields with R c > a/4 the electrostatic correc- 
tions to the bandwidth at the Fermi level become much 
less pronounced and the validity of (|2f| ) is eventually 
restored well before the first magnetic flat-band situa- 
tion. The main reason for this is that now the Hartree 
potential is able to affect the density independently of 
the dispersion at the Fermi level. Consequently, self- 
consistency can be achieved without changing |Ae„ F |. 



We first discuss the weaker modulation without band 
overlaps (Fig. ||). Here, for the higher temperature, the 
density modulation produced by the periodic magnetic 
field is already small without inclusion of the Hartree po- 
tential so that the self-consistent potential is also minute. 
But also for the lower termperature in (b) the situation 
changes around R c — a/4 and the bandwidths become 
close to the noninteracting values, although the ampli- 
tude of the self-consistent potential remains appreciable. 
Around R c — (3/8) a an electric flat-band condition is 
satisfied and the first Fourier component of the Hartree- 
potential has no effect on |Ae„ F |. Therefore here the 
change of the bandwidth at the Fermi level due to the 
Hartree-potential must be small but this is not reflected 
in the amplitudes of the self-consistent densites and po- 
tentials. For still lower fields, away from the electric flat- 
band situation, the Hartree potential yields again notice- 
able corrections to the bandwidths but the noninteract- 
ing curve remains essentially valid. 

For the stronger modulation displayed in Fig. (^), the 
noninteracting bandwidth at the Fermi level is recovered 
as soon as it gets larger than hug. Since the bands at the 
Fermi level then overlap, the density modulation consists 
mainly of higher Fourier components with wave vectors 
r]K, rj > 1. The induced potential therefore belongs al- 
ready to the regime of validity of the linear relation (|2^ ) 
although we still have R c < a/8. Therefore the Hartree 
potential here reduces the density modulation in ampli- 
tude but does not much alter its shape, while the bands 
around the Fermi level remain dominated by the cosine 
form imposed by the periodic magnetic field. 



V. CONCLUSION 

We have calculated the density response of a 2DEG in a 
quantizing magnetic field Bq to a magnetic cosine mod- 
ulation B\ cos(Kx) and compared it with the response 
to an electric cosine modulation. We also included self- 
consistently the induced electrostatic potentials which re- 
duce the density fluctuations. We investigated in detail 
the changes in the energy spectrum brought about by the 
requirement of electrostatic self-consistency. In contrast 
to the case of an electric modulation, where the Hartree 
potential always tends to decrease the width of the Lan- 
dau bands, for a magnetic modulation the Hartree poten- 
tial may either decrease or increase the band dispersion, 
depending on strength and period of the modulation and 
on the filling factor. 

In any case the produced density modulation depends 
crucially on the temperature. If ksT is not small com- 
pared to Hloq the behaviour of the density is quasi- 
classical, i. e. the periodic magnetic field does not lead 
to an appreciable density modulation, while the density 
modulation induced by an electrostatic potential essen- 
tially follows the Thomas-Fermi formula. For temper- 
atures satisfying kgT <ti hujQ (which we assume in the 
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remainder of this section) both types of modulation lead 
to an appreciable non-classical inhomogeneity of the den- 
sity and thus to a non-trivial electrostatic self-consistency 
problem. 

We were mainly interested in the effect of the Hartrec- 
potential on the spectrum around the chemical potential. 
We found that an important parameter is the ratio of 
the cyclotron radius to the length scale a Sp of the den- 
sity variation. If the bands in the vicinity of the Fermi 
level do not overlap we have a Sp ~ a , whereas for over- 
lapping bands the density fluctuations consists of higher 
Fourier components and a 5p is significantly smaller than 
the period a of the modulation. Concerning electrostatic 
effects we can clearly distinguish two regimes by the con- 
ditions R c -c a Sp /A and R c > a Sp /A. This defines for fixed 
total density a distinction between high and lower aver- 
age magnetic fields. The value of Bq around which the 
regime changes depends for a weak modulation only on 
the period while for a sufficiently strong modulation the 
transition takes place when the bands start to overlap. 

For a Sp much larger than the cyclotron radius, i. e for 
high enough Bq, the dispersion of the energy levels is 
changed drastically by the inclusion of electrostatic self- 
consistency. As in the case of a purely electrostatic moi, 
ulation, for which similar screening effects are known,tZl 
the effective Landau bands may be pinned to the Fermi 
level over regions comparable to the period. The nonuni- 
form magnetic field affects the density also in regions 
where the chemical potential lies in a gap between two 
bands. In those regions the density is not constant but 
reproduces the profile of the magnetic field, since the lat- 
ter alters the number of states in the vicinity of each 
center coordinate. 

If the cyclotron radius is not small enough against a Sp , 
Rc > Osp/4, the inclusion of electrostatic self-cosistency 
does not lead to an appreciable change of the dispersion 
of the bands around the Fermi level. The behaviour of 
the latter in the regime where R c ~ a can therefore safely 
be calculated from the noninteracting system. 
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FIG. 1. Density p(x) in units of 1/2itIq for a mag- 
netic modulation of amplitude Bi = 0.1 T in average fields 
Bo = 2.5 — 1.7 T corresponding to filling factors v — 4 — 6. 
The temperature is T = 1 K. The dashed lines are for the 
noninteracting system; the results for v = 4 and v = 6 are 
marked with circles. The lines between these two show the 
successive filling of the n — 2 level when v is increased in 
steps of 1/3. The solid lines display the results for the same 
filling factors after establishing electrostatic self-consistency. 

FIG. 2. Results for a magnetic modulation of ampli- 
tude B x = 0.1 T at (a) v = 5 (Bo = 2.0 T), (b) v = 4 
(B = 2.5 T) and (c) v = 14.3 (B = 0.7 T). The upper 
panel displays the spectra (dashed lines: noninteracting, solid 
lines: self-consistent solution at T = 1 K, dash-dotted lines: 
self-consistent solution at T = 0.1 K; the horizontal straight 
line with dots indicates the position of the chemical potential 
which, to the accuracy of the figure, is the same in all three 
cases) . The lower panel shows the density fluctuation in units 
of 1/27tZq for the self-consistent situations (the solid line is for 
T = 1 K and the dash-dotted line for T = 0.1 K). 



8 



FIG. 3. Bandwidth at the Fermi level and amplitude of 
the induced density fluctuations for the noninteracting sys- 
tems under a weak magnetic modulation Bi — 0.01 T and an 
electric modulation Vi = 0.27 meV (leading to comparable 
maximal bandwidth at the Fermi level, i. e. Vi — A m s) for 
average fields 10 T > B > 0.125 T. The data are plotted 
versus inverse average field indicated on the top x-axis. The 
bottom a;-axis displays the ratio 2R c /a tx B^ 1 . (a) shows the 
width of the band at the Fermi level (solid line for the mag- 
netic, dashed line for the electric modulation; the dash-dotted 
line is Tiuio). The amplitude of the density fluctuation in units 
of l/27rZg is shown in (b) for the magnetic and in (c) for the 
electric modulation. The solid lines are for T = 1 K and the 
dashed lines for T = 0.1 K and the values for even-integer 
filling factors are marked with circles or diamonds, respec- 
tively. In (c) the dash-dotted line displays the prediction of 
Thomas-Fermi-theory |Ap| = 2(Vi/[i)po ■ 

FIG. 4. Bandwidth at the Fermi level and amplitude of 
density modulation for the noninteracting system as in Fig. 
^| but under a stronger magnetic modulation of amplitude 
Bi =0.1 T. In (a) the solid line is the width of the band at 
the Fermi level and the dash-dotted line is Uloq. (b) displays 
the amplitude of the density; the solid line is for T = 1 K and 
the dashed line for T = 0.1 K; the circles and diamonds mark 
the values for even-integer filling factors. 

FIG. 5. Self-consistent results for a magnetic modula- 
tion of amplitude B t = 0.01 T for 10 T > B > 0.18 T 
plotted against 2R c /a. The upper two panels show for (a) 
T = 1 K and (b) for T = 0.1 K the bandwidth |Ae„ F | at the 
Fermi level (dashed lines with diamonds) and the amplitude 
of the self-consistent potential |AV H | (dashed lines with cir- 
cles). The filling factor was increased in steps of 1/3 so that 
each symbol corresponds to a calculated value and the lines 
are only guides to the eye. The dash-dotted line in (a) and 
(b) is Twoo and the solid line displays the noninteracting band- 
width at the Fermi level for comparison. In (c) the amplitude 
of the self-consistent densities are shown; here the solid line 
with circles is for T = 1 K and the dashed line with diamonds 
for T = 0.1 K. 

FIG. 6. Self-consistent results for a magnetic mod- 
ulation of amplitude Bi — 0.1 T and average fields 
10 T > Bq > 0.3 T. The other parameters and the mean- 
ing of the lines is the same as in Fig. n. 
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